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Abstract 



We construct globally regular gravitating solutions, which possess only discrete symme- 
tries. These solutions of Yang-Mills-dilaton theory may be viewed as exact (numerical) solu- 
tions of scalar gravity, by considering the dilaton as a kind of scalar graviton, or as approximate 
solutions of Einstein- Yang-Mills theory. We focus on platonic solutions with cubic symmetry, 
related to a rational map of degree N = 4. We present the first two solutions of the cubic 
N = 4 sequence, and expect this sequence to converge to an extremal Reissner-Nordstrom 
solution with magnetic charge P = 4. 

1 Introduction 

O, 

Localized 'particle-like' solutions play an important role in non-linear field theories. In flat space, 
solutions of particular interest are Skyrmions (in the strong interaction) £Q , electroweak sphalerons 
\ [2], or magnetic monopoles 

In the case of pure Yang-Mills theory, or of pure Einstein gravity, in contrast, there are no 
Q-i, regular 'particle- like' solutions jlj. But when Yang-Mills theory and Einstein gravity are combined, 

static, gobally regular 'particle- like' solutions appear 0. 

The spatial symmetries of static classical solutions can be specified by means of certain rational 
maps of degree TV [B] . Maps of degree N = 1 and N — 2 lead to classical solutions with spherical 
symmetry and axial symmetry, respectively, while for maps with degree N > 3, solutions with only 
discrete symmetries appear. Identifying these with the symmetries of platonic solids or crystals, 
we refer to them as platonic solutions. 

In flat space platonic solutions, such as platonic Skyrmions, platonic monopoles and platonic 
sphalerons, are known [7||H1III]. But does gravity also allow for the existence of regular solutions 
with only discrete symmetries, and does it even allow for platonic solutions, which do not possess 
a flat space limit? 

To answer these questions, we would like to construct platonic generalizations of the famous 
static, spherically symmetric Bartnik-McKinnon (BM) solutions 0, representing globally regular- 
solutions of Einstein- Yang-Mills (EYM) theory with action 

s = / (f " ^ V")) V~gd* x , (i) 

where R denotes the scalar curvature, F^ v = d^A v — d v A^ + i \A^,A V \ the SU(2) field strength 
tensor, A^ = A°r a /2 the gauge potential, and Newton's constant G = 1/4tt. 

For static solutions, the 4-dimcnsional metric (x^ = t,x l ) can be reduced according to 

ds 2 = -e 20 di 2 + e-^dsj , ds 2 = h ij dx i dx j , (2) 
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yielding the 3-dimcnsional action 



S = y ~ l d ^ d ^ ~ \^ VM 3 x , (3) 

with A^dx^ = Aidx 1 . 

While we expect a 3-dimensional metric hij — hi(x)Sij with distinct metric functions hi for 
exact platonic solutions, we here consider only approximate platonic EYM solutions, obtained 
with a euclidian 3-dimensional metric. For static configurations, in this approximation the action 
then agrees with the action of Yang-Mills-dilaton (YMD) theory, 

S = J {-^-^^n}^ . (4) 

where 4> denotes the dilaton field, and the solutions therefore represent platonic solutions of YMD 
theory. Indeed, the dilaton may be considered as a kind of scalar graviton, possessing also a 
universal coupling to matter. 

Consequently, solutions of EYM and YMD theory are expected to possess very similar prop- 
erties. This is indeed the case for spherically symmetric and axially symmetric solutions 1101 
1111 IT^j . The study of platonic YMD solutions is therefore expected to yield important informa- 
tion about platonic EYM solutions, without having to solve the highly complicated set of coupled 
Einstein- Yang-Mills equations. 

Spherically and axially symmetric EYM and YMD solutions form sequences, labelled by two 
discrete numbers, representing the degree of the map N, and the node number k of the gauge field 
functions (HlEB^DEl- We expect analogous sequences of platonic solutions, i.e., for a given map 
of degree N, there should be a fundamental platonic solution and an infinite sequence of excited 
platonic solutions with the same symmetries. First indications of the existence of these solutions 
have been seen recently in Yang-Mills-Higgs-dilaton (YMHD) theory 

Here we construct such platonic YMD solutions. We consider the Ansatze for these solutions 
in section 2, we discuss their properties and present the numerical results in section 3, and we give 
the conclusions and an outlook in section 4. 



2 Ansatze for Platonic Solutions 

To obtain YMD solutions with discrete symmetry we make use of rational maps, i.e. holomorphic 
functions from S 2 i— > S 2 [jj]. Treating each S 2 as a Riemann sphere, the first having coordinate £, 
a rational map of degree N is a function R : S 2 i— > S 2 where 

jko = § , 0» 

and p and q are polynomials of degree at most N, where at least one of p and q must have degree 
precisely N, and p and q must have no common factors jS] . 

We recall that via stereographic projection, the complex coordinate £ on a sphere can be 
identified with conventional polar coordinates by £ = tan(6*/2)e lv 0. Thus the point £ corresponds 
to the unit vector 

n f ~(23}(Q,23(0,l-|a. (6) 
and the value of the rational map R(£) is associated with the unit vector 

%. = TT ^(25R( J R),23(i?),l-|i?| 2 ) . (7) 

We here consider platonic YMD solutions obtained from the map R4, 

*(0 - e4 + 2 f^ + 1 • (8) 
V ' £ 4 -2 v / 3ie 2 + l 
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This rational map leads to the unit vector (JZJ) 

/ (r 2 - z 2 ) 2 - 2z 2 r 2 + 2x 2 y 2 VS(r 2 + z 2 )(x 2 - y 2 ) Ay/Zrxyz\ 
nR * = \ A/4 ' M '— ~ J ' (9) 

with A/4 = 2(x 4 + x 2 y 2 + x 2 z 2 + y 4 + y 2 z 2 + z A ). 

For static solutions with Aq — 0, Ai — Ai(x), i = 1,2,3, 4> = 4>{x), the general set of field 
equations then involves the dilaton function 4>(x), and nine functions for the gauge field, A®(x). 
The solutions are characterized by the rational map of degree N, and, for a given map, numbered 
with increasing excitation by the integer k, where the fundamental solution has k = 1. 

To obtain solutions with finite energy, we demand that the gauge field of the fundamental 
solution (and for odd k solutions in general) tends to a pure gauge at infinity, 

Ai = {(diU^Ul , (10) 

with Uoo — ir R , where 

tr := n R ■ f . (11) 

In contrast, the gauge field of the first excited solution (and for even k solutions in general) should 
vanish at infinity. 

We also require that the dilaton field vanishes at infinity, 0(oo) = 0, since any finite value of 
the dilaton field at infinity can always be transformed to zero via <p — > <f> — 4>(oo), when the radial 
coordinate is transformed according to r — > re~^°°\ 

The boundary conditions for platonic solutions are, in principle, only needed at infinity. Speci- 
fying the unit vector n R = n Ri thus completely determines the boundary conditions for the gauge 
field of the fundamental N — 4 platonic solution with cubic symmetry p0(l . The gauge field then 
exhibits at infinity the angle-dependence of the rational map, since Uoo = irn 4 . In contrast, for 
the first excited N — 4 platonic solution the discrete symmetry does not enter via the boundary 
conditions at infinity, since the gauge field vanishes there. 

Parametrizing the gauge potential Af in terms of polar components, Ag, and A", we 
obtain a smooth numerical solution for the first excited platonic N = 4 YMD solution. For 
the corresponding fundamental YMD solution, however, this parametrization does not lead to an 
everywhere smooth numerical solution (with the same numerical method |14|). The polar gauge 
field functions of the fundamental solution exhibit a discontinuity at infinity, which we relate to 
the fact that the polar gauge field functions of the fundamental solution are angle-dependent at 
infinity, whereas those of the first excited solution approach zero at infinity. 

Assuming that the observed discontinuity is a numerical effect, we expect to be able to avoid 
it by employing a better suited parametrization for the gauge field. In particular, the angular 
dependence of the gauge field at infinity should be contained explicitly in the new parametrization 
for the fundamental solution [15]. 

To obtain such a parametrization, we note, that at infinity Ag and A v approach 

i {dgU^Ul = , and * (d^U^Ul = -sinfl d e r R , (12) 

respectively, since 

n R x d e n R — d v n R / sm.6 . (13) 

We therefore choose t r , dgr R , and d v T R as a basis of the Lie algebra of SU(2). This is adequate, 
as long as dg n R and d v ri R do not vanish. Since for the rational map -R4, the vectors dg n Ri and 
dip n Rl vanish along the cartesian coordinate axes, problems may arise in the numerical proceduce, 
unless boundary conditions are specified here. Thus we demand that the gauge field vanishes along 
the coordinate axes. 

Accordingly, we employ the following parametrization for the gauge field of the fundamental 
solution, 

. f l d v T R 2 dgT R 3 

A e = Jg -r-^ + Jg m „ , + Jg T R , (14) 
sm0 \ogn R \ 
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xz- and yz-plane (93 = or 7r/2) 


ip = 7r/4 


1 


# = , # = , 9^ = 0, 
^/e=0 , /| = , .f e 3 = , 
d<pfi = , f*=0 , d v fl = , 


# = , d v tf = o, # = o, 

d v fe=0 , fe=0 , f$ = , 
d v ti = , /*=0 , 9 v / 3 = , 


2 


A$ = , >4£ = , d^A* = , 
i4j = , A* = , 0^ = , 
9^ = , 9^ = ; ^ = , 


^ = , 9 v Aji = , ^ = , 
4J = , 9 V A* = , Ag = , 
9 V A*=0 , A*=0 , 9^ = , 




z-axis (9 = 0) 


a;y-plane (0 = 7r/2) 


1 


# = , # = , d e f ? = , 
defi = , f 2 e =Q , / e 3 = , 

fc/£ = o , /| = o , /3 = o , 


/ r A = , f? = , d e f? = , 
doft=0 , f$ = , dgf$ = , 
defi = , e = , /| = , 


2 


^ = , Aj? = , 9 e ^ = , 

= , A 2 e = , ^| = 0, 
AJ, = , A*=0 , A 3 =0 , 


A* = , Af. = , 9 e A 3 = , 
fl i4j = O , d B A 2 e =0 , A 3 e = , 
Al = Q , Al = , 9^ = . 



Table 1: The boundary conditions in the xz—, yz — , and sy-plane, on the z-axis, and for tp — n/4 are 
given for the gauge field components of the fundamental (k — 1) and first excited (k = 2) platonic solution. 



A v = 4 (- sin 6 dgr R ) + % + / 3 r R , (15) 

v \o v n R \ 

and 

A r = frT x + f?T y +f?T z . (16) 

The functions /" now have no angular dependence at infinity. In particular, fg and ft approach 
unity at infinity, while the other seven functions vanish there. 

Subject to the boundary conditions and to the gauge condition 

d l A = , (17) 

the set of field equations can be solved numerically. The discrete symmetries of the platonic 
solution allow for a restriction of the region of numerical integration. Specifying the reflection 
symmetries of the fields with respect to the xy-, xz- and yz-plane, it is sufficient to solve the field 
equations for x > and y > and z > only |13j . The integration region can even be restricted 
to x > y > and z > 0, with the reflection symmetries of the fields at ip — it /A . 

The boundary conditions at infinity must now be supplemented by conditions at the other 
boundaries of the integration region. These boundaries consist of the origin, parts of the xy-, 
xz- and yz-plane (or the plane tp = 7r/4), and the positive z-axis. (Together with the boundary 
at infinity these boundaries are implemented numerically as the boundaries of a cube |14|.) The 
gauge field vanishes at the origin, and the dilaton field satisfies d r (f> = 0. 

The boundary conditions in the a;z-plane (tp — 0), the yz-plane {ip = tt/2), in the plane tp — 7r/4, 
on the z-axis (6 — 0) and in the xy-plane (6 = tt/2) follow from the reflection symmetries of the 
fields. In particular, we suppose that the rational map R4 determines the reflection symmetries 
of the fields. Thus we impose, that for the platonic YMD solutions the gauge field Ai possess the 
same reflection symmetries as i[9jTR 4 , trJ. For the gauge field functions Af and the proper 
set of boundary conditions is given in Table 1. For the dilaton field the normal derivative must 
vanish in the xy-, xz- and yz-plane, in the plane tp = 7r/4, and on the z-axis (dg<p = 0). 

An adequate initial guess is mandatory for the numerical construction of platonic solutions. 
The initial guess is obtained by parametrizing the gauge field according to 

1 - fir) 

Ai= ^(d t r Ri )r Ri , (18) 

and the dilaton field by 4>(r). This simple ansatz with two radial functions f(r) and <j>(r) is inserted 
into the action, the angular dependence is integrated, and a set of variational equations is obtained, 
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.829 1.385 1.870 2.319 
.971 1.796 2.527 3.197 
.995 1.935 2.805 3.621 
.999 1.980 2.921 3.824 


.804 1.336 1.800 2.231 
.966 1.773 2.482 3.137 
.994 1.927 2.785 3.588 
.999 1.977 2.913 3.808 
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Table 2 

The energy E of EYM solutions and YMD solutions with spherical symmetry (N = 1) and axial 
symmetry (N = 2 — 4) and node numbers k = 1 — 4. For comparison, the energy Eoo of the 
limiting solutions of the sequences is also shown. Note, that E/N decreases with N for fixed finite 
k. 



similar to the set of spherically symmetric (N = 1) equations. These ordinary differential equations 
are solved numerically, subject to the set of boundary conditions of the spherical (N = 1) solutions 
(in particular /(0) = 1, /(oo) = —1 for the fundamental solution and /(0) = 1, /(oo) = 1 for the 
first excited solution). The initial guess is then given by (|18f) with the numerically determined 
functions f(r) and <f>(r). 



3 Properties of Platonic Solutions 

The numerical solutions are constructed with help of the software package FIDISOL ^1] based 
on the Newton-Raphson algorithm. The solutions are obtained in spherical coordinates r, 6, (p. 
To map the infinite range of the radial variable r to the finite interval [0, 1] we introduce the 
compactified variable r = r/(l + r). Typical grids contain 70 x 25 x 25 points. The estimated 
relative errors are on the order of 1%, or better. 

For static finite energy solutions (with Aq = 0) the energy is given by 

E = J Qdi W + ~e 2 *Tr (FyF«)) d 3 x . (19) 
It is related to the dilaton charge D, 



D = TZ I V<t>-dS, (20) 
Is-, 



4tt 



by nm 

The energy density e(x) is defined by 



D = E. (21) 



E = J e{x)d 3 x . (22) 



Let us begin the discussion of the properties of the platonic YMD solutions by briefly recalling 
the properties of the spherically and axially symmetric EYM and YMD solutions [51 1101 ITT1 112| . 
In both cases for an appropriate map of degree N, an infinite sequence of solutions exists, labelled 
by the node number k of the gauge field function(s). 

The energy E of the solutions of the sequences corresponding to maps with N = 1 — 4 is shown 
in Table 2 [111 I12j . and exhibits close agreement between the sets of EYM and YMD solutions. 

In the limit k — > oo, the sequences of solutions corresponding to a certain map of degree N tend 
to limiting solutions. For the EYM sequence, the limiting solution corresponds to the exterior of 
an extremal Reissner-Nordstrom (RN) solution with mass M = N, and magnetic charge P — N. 
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In isotropic coordinates, the horizon of the extremal RN is located at rjj = 0, and its metric reads 



ds 2 = - (j^^j dt2 + (^7^) ( dr2 + 7-2 ( d ° 2 + sin2 ed v 2 )) • ( 23 ) 

This metric is of the type of the metric Eq. J5J) with 



N + r 



(24) 



and a euclidian 3-dimensional metric. Not surprisingly, for the YMD sequence the limiting solution 
has dilaton field [HI E] 

^ = -lnfl+^j , (25) 

clearly representing the same extremal RN solution. Thus with increasing excitation EYM and 
YMD solutions differ less and less, and converge to the same limiting solution. Considering a 
YMD solution as an approximation to an EYM solution, the approximation thus improves with 
increasing k for fixed N. 

Let us now turn to the platonic YMD solutions. In Fig. 1 we present surfaces of constant total 
energy density etot and of constant energy density of the gauge field £f = ^e 2< ^Tr (F^F^) for the 
fundamental and first excited cubic solutions. The energy density clearly reflects the symmetries 
of a cube. As before j^lEI we conclude that the shape of the energy density of platonic solutions is 
determined primarily by the rational map, since we observe the same shape for the YMD solutions 
as for Skyrmions, sphalerons, and monopoles, based on the same map R4 0E1IE]- The energy 
density of the excited solution exhibits a cubic shape both for large and for small constant values, 
revealing a small interior cube within a larger exterior cube. 

The energy of the fundamental cubic solution is Ek=i — 2.203, while the energy of the first 
excited cubic solution is E^=2 = 3.11, with estimated errors on the order of 10 -3 and 1CP 2 , 
respectively. Comparison of these values with those of the axially symmetric N = 4 solutions 
(Table 2) reveals, that the platonic solutions have lower energy. 1 For the cubic EYM solutions 
we conjecture that they possess a slightly higher mass, in accordance with the known rotationally 
symmetric solutions. The observed values Ek=i and E^=2 are in agreement with the conjecture, 
that the energy of the excited cubic solutions tends to the limiting value E^ = 4 for k — > 00. 

The metric function —goo = e 2 * is exhibited in Fig. 2 for the fundamental and the first excited 
solution. Clearly, the metric exhibits the same platonic symmetry as the energy density. 

In Fig. 3 we further address the conjecture, that platonic solutions associated with maps of 
degree N form infinite sequences, tending to extremal RN solutions with magnetic charge P = N. 
The cubic solutions should therefore tend to a limiting extremal RN solution with charge P = 4. 

We first consider convergence of the dilaton function 4>k(x) to the limiting function (j)oc(x), 
Eq. 1)25(1 . Therefore we extract a 'local magnetic charge' Pk(x) for the fundamental (k = 1) and 
first excited (k — 2) solution, by defining 

MS ) = -In ( 1 + ■ (26) 



Exhibiting Pk{x) as a function of r for three fixed angles, we observe indeed that the functions 
Pk (x ) approach the constant P^ — 4 with increasing k in an increasing interval. 

We then extract a 'local magnetic charge' Pk{x) from the metric component goo of a non- 
extremal RN solution of mass E in isotropic coordinates [H], employing Eq. J2J) 

- ^ G + e-S + S + f l + e-^) ) • (27) 



1 Note that the difference in energy for the first excited solution is the same order of magnitude as the estimated 
error though. 



() 



k=1 jGlof 0.4 k=1,e ^0.002 




Fig. 1 Surfaces of constant total energy density e t ot (left column) and of constant energy density of the 
gauge field ef (right column) are shown for the fundamental cubic YMD solution (upper row) and for the 
first excitation (lower row). 



Again we observe, that the charge P2 (x ) obtained from the first excited solution is closer to the 
constant value — 4 of the limiting extremal RN solution in a wider range than the charge 
P\ (x ) of the fundamental solution, indicating that with increasing k the platonic solutions indeed 
tend to the spherically symmetric extremal RN solution with P = 4. 

4 Conclusions and Outlook 

Motivated by the aim to demonstrate the existence of non-perturbative gravitating solutions with 
only discrete symmetries, we have constructed platonic solutions in YMD theory. These solutions 
may be viewed as exact (numerical) solutions of scalar gravity, by considering the dilaton as a kind 
of scalar graviton, or as approximate solutions of Einstein- Yang-Mills theory, when the spatial part 
of the metric is taken to be euclidian. 

We have focussed on platonic solutions with cubic symmetry, related to a certain rational map 
of degree N = 4. By constructing numerically the fundamental cubic solution and its first excita- 
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k= 1 ; -g 00 =0 .07979 k=2, -g oo =0 .019125 




Fig. 2 Surfaces of constant metric function —goo = e 2 ^ for the fundamental cubic YMD solution (left 
column) and the first excitation (right column). 



tion, we have here presented the first two solutions of the cubic N = 4 sequence. We conjecture, 
that the sequence converges to the extremal Reissner-Nordstrom solution with magnetic charge 
P = 4. 

By including the dilaton to approximate the effects of gravity, we have obtained evidence 
for the existence of platonic EYM solutions, while avoiding the complexity of the full set of EYM 
equations in the absence of rotational symmetry. Construction of exact (numerical) platonic EYM 
solutions, however, remains a true challenge. 

The existence of gravitating regular solutions involving non-Abelian fields, is related to the 
existence of black holes with non-Abelian hair JT^l- Thus the existence of gravitating regular 
solutions with platonic symmetries strongly indicates the existence of a completely new type of 
black holes: static black holes which possess only discrete symmetries |17| . 
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Fig. 3 The 'local magnetic charge' Pk(x) extracted from Eq. (left) and Pk{x) extracted from Eq. ^21\ 
(right) for the fundamental (k — 1) and first excited (k — 2) cubic solution for three representative angles. 
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